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1.  INTRODUCTION 


The  problem  of  testing  for  the  dimensionality  of  the  regression  matrix  under 
multivariate  regression  model  has  received  considerable  attention  in  the  litera¬ 
ture.  This  problem  arises  in  the  areas  of  pattern  recognition,  signal  processing,, 
and  functional  and  structural  relations.  For  a  discussion  on  applications  in 
functional  relations,  the  reader  is  referred  to  Anderson (1 984) . 

Fisher  (1938)  considered  the  problem  of  testing  for  the  number  of  significant 

I 

discriminant  functions  and  it  is  a  special  case  of  the  problem  of  testing  for  j 

the  rank  of  the  regression  matrix.  Tintner  (1945)  derived  the  likelihood  ratio 
test  (LRT)  statistic  for  the  rank  of  the  above  matrix  when  the  covariance  matrix 
is  known.  Anderson  (1951)  derived  the (LRT)  statistic  for  the  rank  of  the  regres¬ 
sion  matrix  when  the  covariance  matrix  is  known.  Fujikoshi  (1974)  derived  the 
LRT  procedure  for  the  rank  of  the  regression  matrix  under  growth  curve  model. 

Recently,  Rao  (1985)  considered  the  LRT  procedure  for  the  rank  under  a  general 

1 

model,  incorporating  the  multivariate  regression  model  and  the  two-way  classifica¬ 
tion  with  interaction  and  with  one  observation  per  cell.  The  above  work  was  done 
when  the  underlying  distribution  is  multivariate  normal.  The  object  of  this  paper 

is  to  discuss  various  procedures  for  testing  for  the  dimensionality  of  the  regression 

. 

matrices  and  derive  asymptotic  distributions  of  the  test  statistics  when  the 
underlying  distribution  is  real  or  complex  elliptically  symmetric  distribution. 

In  Section  2  of  this  paper,  we  give  some  preliminaries  and  state  the  main  j 

problems  that  are  considered.  The  LRT  procedures  for  the  dimensionality  of  the  j 

regression  matrices  are  derived  in  Sections  3  and  4  for  the  cases  of  the  real 
and  complex  elliptically  symmetric  distributions  respectively.  Asymptotic 
distributions  of  the  above  test  statistics  are  derived  in  Section  5  when  the 


joint  distribution  of  the  observations  is  elliptically  symmetric.  Multivariate 
normal  and  multivariate  t  distributions  are  special  cases  of  the  elliptically 
symmetric  distributions.  In  Section  6,  we  derive  the  asymptotic  distribution 
of  the  LRT  statistic  for  the  rank  of  the  regression  matrix  when  the  observations 
are  distributed  indepently  as  elliptically  symmetric.  The  assumptions  made  about 
the  underlying  distributions  is  Sections  5  and  6  are  equivalent  only  in  the  case 


of  multivariate  normal. 


Accession  For 


2.  STATEMENT  OF  PROBLEMS  AND  PRELIMINARIES 
Consider  the  model 


X  -  AO  +  E 


(2.1) 


where  the  error  matrix  E:  n*p  is  distributed  as  an  elliptically  symmetric  dis¬ 
tribution  with  density 


f (E)  =  — hCtrSfVE) 

Mn/ 


(2.2) 


where  h(x)  is  strictly  decreasing  and  differentiable  function  of  x.  Also, 
A:  n*k  is  known  and  of  rank  k<  n,  and  8:  k*p  is  unknown.  Now,  let 


A  *  CO 


(2.3) 


where  C:  uxk  is  known  and  of  rank  u.  Let  denote  the  hypothesis 
that  the  rank  of  A  is  r  and  denote  the  hypothesis  that  Ae  P^.  Here  Ae  P^. 
denotes  that  the  rows  of  A  lie  in  a  r-dimensional  plane  in  p-dimensional  space. 
Now,  let  n^(a)  denote  the  set  of  n*p  matrices  of  the  form 


M-  (GF+  ab  '  )D 


C2.4) 


where  | G *G |  0 ,  FF**  1^  and  D:  pxp  is  any  positive  definite  matrix  and  b  is  any 

p*l  vector.  Then  H^  denotes  the  hypothesis  that  Ae  ifr(0)  whereas  ^  denotes  the 


hypothesis  that  Ae  1^(1)  where  1'  ■  (1,...,1). 
Next,  consider  the  model 


Z  -  AS  +  N 


(2.5) 


where  8:  kxp  is  an  unknown  complex  matrix,  N:  n*p  is  distributed  as  the  complex 
elliptically  symmetric  distribution  with  density  given  by 


& 

- -  hC2tr  E-1N'N)  . 

ur 


(2.6) 


Here,  we  note  that  the  complex  elliptical  distribution  was  introduced  by  Krishnaiah 
and  Lin  (1984).  In  this  paper,  we  consider  the  problems  of  testing  the  hypotheses 
Hj  and  H2  when  the  «sad«rlying  distributions  are  real  and  complex  elliptically 
symmetric.  We  need  the  following  lemmas  in  the  sequel. 

Lemma  2.1.  If  a  is  a  u*l  complex  vector,  L:  nxp  is  a  complex  matirx,  Q:  pxp  is 
Hermitian,  positive  definite  matrix  such  that  the  rank  of 

S  -  Q_1/2L'(I  -aCa,a)“1a,)LQ_1/2 


is  r  or  more,  then  the  eigenvalues  of 


-1/2  -  -1/2 

S (M)  =  Q  '  (L-M)‘(L-M)Q 


are  minimized  simultaneously  with  respect  to  Me  ir^a)  if  and  only  if 


M=  [a(a'a)-1a'LQ-1/2+(I  -a(.a*a)  'a' )LQ~1/2V’VjQ  1/2 


n  _  -  _ 


where  the  rows  of  consist  of  normalized  eigenvectors  corresponding  to  the  first 

A 

r  largest  eigenvalues  of  S.  The  minimum  values  of  ch^(S(M))  are  given  by  4>r+i 
where  ch^A)  denotes  i-th  largest  eigenvalue  of  A,  .  .^4^  are  the  eigenvalues  of 

A 

S  and  j  *  0  for  j^_p. 

When  a,  L  and  Q  are  real,  the  above  lemma  was  proved  by  Fujikoshi  (1974).  The 
proof  of  Lemma  2.1  follows  along  the  same  lines  as  in  Fujikoshi  (1974).  From 
Lemma  2.1,  the  following  lemma  follows  immediately: 

Lemma  2.2.  Let  f(S)  be  a  function  of  a  pxp  Hermitian  matrix  S  such  that 


f(S)  *»  g(ch,  (3)  ... .  ,ch_(S)) 


5 


31 


% 

is 


a 


and  g(*)  is  strictly  decreasing  in  each  argument.  Then 

max  f(S(M))-g(<i  . A.  0,...,0) 

Me*  (a)  r+1  p 

r  _ 

where  <frj  -  ch^ (S)  for  j  =  l,2,...,p. 

The  following  lemma  was  proved  in  Bai  (1984). 

Lemma  2.3.  Suppose 

f  (7\  -  aCn>71C4.a(n>7K^14.  4.  (n) 

fn(Z)  -  aK  Z  +aK_xZ  +.,.+a0 

f(Z)  ■akZk+ak.1Zk-1+...+ao 

where  K>_k.  Also,  let  fn<Z)->f(Z)  as  where  a^O,  n-1,2,..,  and  a^O. 

In  addition,  let  Z^,...,Z^  denote  the  roots  of  f  (Z)  .  Then,  we  can  suitably  arrange 
the  roots  of  fQ  as  Z^n) , . . . ,Z^n) , . . . ,Z^n)  such  that 

Zi(n)  +  Zi  for  i<^k 
|zj  -  00  for  i>  k 

as  n ->°°. 


'  "■  ■  *  ■  *v_v>*.yy; -.V.V. 


3.  LR  TESTS  FOR  THE  DIMENSIONALITY  OF  THE  REGRESSION  MATRIX 


In  this  section,  we  derive  the  LR  test  for  the  dimensionality  of  the  regression 

matrix  under  the  model  (2.1)  for  the  cases  when  Z  is  known  and  unknown  and  the 

underlying  distribution  is  elliptically  symmetric.  For  the  sake  of  simplicity,  we 

first  reduce  the  model  to  a  canonical  form.  It  is  known  that  nonsingular  matrices 

Ta,  Tg  and  orthogonal  matrices  TA»  exists  such  that 
k*k  u*u  n*n  k*k 


A  -  r  YVVr. 

n*k  ^\0  /  ^ 

T  (I  0)r  . 

C  u  C 


(3.1) 

(3.2) 


We  make  the  following  orthogonal  transformation 


Y  -  E  - 


rc  ° 


n-k 


r!  (X  -  A0)  . 
A 


From  assumption  (2.2)  we  get  the  density  function  of  Y  as 

n 

|Z|  2  h(tr!_1(Y-  E)  '  (Y-  E)  )  . 


(3.3) 


(3.4) 


If  we  partition  Y  and 


Y  = 


-  as 


f 

Y1 

u 

'1 

Y2 

k-u 

z2 

Y3J 

n-k 

"3 

P 

P 

it  is  easy  to  see  that 


V- 


yi  ■  (iuo)rc(ik  o)r;x 


Y2  -  <°  Wrc(Ik  0)rAX 


Y3  ■  (0  W? 


H1  ■  (In  °>WKA6  -  Va 


(3.5) 


“2  ‘  <°  Wc(Ik0)rlAe=  t°  WVa9 
S3  *  <0  Wi*8  ■  (°  W  (ok)V  *  °- 

Under  canonical  form  the  hypotheses  and  H 2  are  equivalent  to  E^e  it  (0)  and 
=  !«  irr(a)  respectively  where  a  =  Tc  1.  Now,  let  ^  and  H2  denote  the  alternative 
hypotheses  E^e  7^,(0)  and  E^  e  TTr,(a)  respectively  for  some  r'  >  r.  Also,  let 

M  -  CCA’ArV 


Eq  *  (A1 A)_1A'X 
Sh(5,M)  =  (C  5) ,M~1(C5) 


(3.6) 


S  (5,M)  «  (CS),{m“1-m”11(1*m“11)”11»m“1}CE 


S  =  X'(I-A(A,A)~1A,)X. 


3 . 1  LRT  Statistics  when  Z  is  Known 


When  Z  is  known,  the  LRT  statistic  for  testing  against  is  given  by 

^  1 

max  sup  Izl"  7  h(trZ(Y-E)'(Y-E)) 


7  h(trZ*-L(Y-E)  '(Y-E)) 


max  h(trE  (Y 


)+trE  YIYJ 


h(trE~  i^Y-j) 


Since  h(x)  is  a  decreasing  function  of  x,  Y^£  ^Y^>_  0  and  (Y  -E^)£  3(Y^-5)'>  0, 


we  obtain,  from  Lemma  2.1, 


h(<l>r+l+---+4>s+trE":LY3Y3) 


h(trE_1Y^Y3) 


(3.7) 


where  Y^  =  X'd^-  A(A' A)_1A' )X,  >0  are  positive  eigenvalues  of 

-1/2  -1/2  -1 

1  YjY^E  =  Sh^“0’^E  and  S  =  m^n(u>P)"  When  the  underlying  distribution 
is  multivariate  normal,  we  obtain 


-2  log  Tx  -  ♦r+1+. .  .ns 


(3.8) 


The  LRT  statistic  for  testing  against  H2  is  given  by 

max  hCtrE"1(Y1-H1),(Y-51)+tr£"1YlY,) 

=.€*  (a)  11  11  33 

l  r 


h(trl“1Y^Y3) 


h(<l'r+i+- +  trE_1Y^Y3) 


h(trE_1Y^Y3) 


(3.9) 


where  0  are  the  positive  eigenvalues  of  S^(5q,M)E  and  s  =  min(u-l,p) 

When  the  underlying  distribution  is  multivariate  normal,  we  obtain 


-2  log  T2  *  W**4^  • 


(3.10) 


3 . 2  LRT  Statistics  when  E  is  Unknown 

* 

When  E  is  unknown,  the  LRT  statistic  for  testing  H,  against  H  is  given  by 


max  max  sup  |z|  h(trZ  (Y-E) '  (Y— IZ ) ) 

SlS*r<0>  E>0  ~2 

h - - - - - 

max  sup  |z  |  ^  h(trZ  i(Y-E) ' (Y-E)) 

l> o  e2 

n_ 

max  max  |Z|  2  h(trZ~1(Y  -E. ) ' (Y  -S)+Y!Y  )) 

-  __  /A\  r i-v  1  1  i.  J  -5 


=€*(0)  Z>0 


max  |z|  2  h(trZ  ^YiY  ) 

Z>0  J 

_  S£ 

max  A  2  (Y1-5- ) ' (Y1-E1 )+Y'Y 

,ns  max(h)  11  11  3  3 

-jtr  co) 

np  n  _ 

A"  2,.  JyIyJ"  2  h(r— 2 - ) 

ma«h)  3  3  Xm«<h) 


max(h) 


(3.11) 


Equation  (3.11)  follows  from  Anderson  and  Fang  (1982)  who  derived  the  LRT  procedure 
for  0=0  when  the  underlying  distribution  is  real  elliptically  symmetric.  Now, 
using  Lemma  2.1,  we  obtain 

.1  .i-f 

T3  Z«VHl)'0rrSl))(Y3Y3)  2 

“le7rr(5) 


=  ((l+dr+1)(l+dr+2)....(l+ds)) 


(3.12) 


where  s  =  min(u,p)  and  d^>_.  .  . >dg  >  0  are  the  positive  eigenvalues  of  S^(Eq,M)S 

* 

Similarly,  the  LRT  statistic  for  testing  H2  against  H2  is 


.  i  .  i  - 1 

T  -  max  |I-(Y'Y  )  2 ( (Y  -3^  '  (Yj-^)) (YjYj)  2 


(3  13) 


[(1+tr+l)‘*,(1+il£)1 


(3.14) 


where  s  =  min(u-l,p)  and  >  0  are  the  positive  eigenvalues  of  Sf(EQ,M)S 

When  the  underlying  distribution  is  multivariate  normal,  the  test  statistics  T 
and  were  derived  by  Fujikoshi  (1974)  and  T2  was  derived  by  Rao  (1965). 


10 


4.  LRT  STATISTICS  FOR  THE  DIMENSIONALITY  OF 
REGRESSION  MATRIX  IN  COMPLEX  ELLIPTICAL  CASE 

Consider  the  complex  multivariate  regression  model  (2.5)  where  N  is  distri¬ 
buted  as  (2.6).  Also,  let 


A0  =  C8  (4.1) 

where  C:  u*k  is  known  and  is  of  rank  u.  In  addition,  let  denote  the 

* 

hypotheses  that  the  rank  of  A^  is  r  whereas  denote  the  alternative  hypothesis 

that  the  rank  of  is  greater  than  r  where  r  <  s  =  min(u,p) .  Also,  let 

* 

denote  the  hypothesis  that  A^  e  and  let  H^q  denote  the  alternative  hypothesis 

that  A.  e  P  ,  for  some  r*  >  r  and  r  <s  =  min(u-l,p).  Here  AeP  means  that  the 

Or  ,r  r 

rows  of  Aq  actually  lie  in  a  r  dimensional  complex  plane.  The  hypotheses 

★ 

and  H20  are  respectively  equivalent  to  A  e  tt^(0)  and  A  e  11^(1).  Also,  let 

* 

denote  the  alternative  hypothesis  that  A  e  r^(0)  for  some  r  >  r  and 

denote  the  alternative  hypothesis  that  A  €  tt^Q)  for  some  r*  >  r.  We  now  reduce  the 

model  in  canonical  form  as  in  the  real  case. 

* 

The  problem  of  testing  H^^  against  H^q  is  equivalent  to  testing  the  hypothesis 
e  xr^CO)  against  the  alternative  e  ir^,  (0)  for  some  r*  >  r  and  r<s  =  min(u,p) 

k 

in  the  canonical  form.  Similarly,  the  problem  of  testing  H  against  H  is 

ZU  Zw 

equivalent  to  testing  the  hypothesis  E^exr^(a)  against  the  alternative 

£  n  ,  (a)  for  some  r'  >  r  and  a  =  T  ^1  and  r  <s  =min(u-l,p).  When  E  is  known, 

1  r  w  «  v  „ 

k 

let  denote  the  LRT  statistic  for  testing  H  against  H  and  let  T,  denote 
j  1U  1U  o 

k 

the  LRT  statistic  for  testing  H2q  against  H^.  Then,  using  Lemma  2.1,  we 
obtain  the  following: 

h (2$  .+. .  .+2*  +2  tr  E_1y;Yo 

T-  =  - £ti - ^  (4.2) 

fi(2  trl  Hf'Y  ) 


W".  "  •  "■  t.  *1  ..  .  .  jr*  ■  .  .  -  .  -  .  ■  ..J'  .  ' 

•'  1*.  n  M  V./ mj?  mjH  m-1? .  f  ~  -V-’'  JL. 


*  *•_ 


m“+  »>'»  I."* 


A  , 


hC2<J<  ,+...+2*  +2  trZ  lYlYj 

- - Eti - _I - U—  (4. 

6  h(2  trZ-iYjY3) 

*  _1 

where  4>  >•••>$  are  the  nonzero  eigenvalues  of  S,  (9 ,M)S.  (I  ,M)  ,  and 

i —  -  s  h  Jc  p  1 —  — 

are  the  nonzero  eigenvalues  of  Sr(0,M)Sn(I  ,M)  \  Here 

i  l  p 

S,  C6,M)  =  (C9)  'M_1CcS) 
n 

W2)  =  z 

M  =  CCA'A}-1^  (4. 

Sf(0,M)  =  CC0),[H"1-M“11C1,M_11)"11,M"1](C0). 

When  the  underlying  distribution  is  complex  multivariate  normal,  we  have 


exp{-(4>r+1+.  ..+$)} 


exp{-  0))^+.  •  .■H'g) }  • 


When  E  is  unknown,  we  denote  the  LRT  statistic  for  and  H^q  against  H 
and  H2Q  by  and  respectively.  Then 


t7-(a^1)awrt2)...aMs))- 


Ts  -  f  u+tj+j)  d+tr+2) . . .  a+«5>  > 

*  -1 

where  d- >..•><!  are  the  nonzero  eigenvalues  of  (0 ,M) {S„  (I ,S) }  and 
l —  —  s  n  x, 

A  -1  • 

are  the  nonzero  eigenvalues  of  (9  ,M) {S^ (I,S) } 


mm, 


5.  ASYMPTOTIC  DISTRIBUTIONS  OF  LRT  TEST 
STATISTICS  FOR  THE  DIMENSIONALITY  OF  REGRESSION  MATRIX 


We  know  that 


Y~  |E|"n/2h Ctrl"1  CY^)’ 0X^2)+ CY2-H2)'CY2-E)+Y'Y3) 


and  d^...>d  >0  are  the  positive  roots  of 


0  =  lWM)-dS|  =  lYiYi  *  dY3Y3l*  0 

Since  H^eir^CO),  i.e.,  rk(E^)*r>0,  we  know  that 

lE^-nAZl-0  (i 

will  have  zero  roots  with  multiplicity  p-r  and  r  nonzero  roots.  We  arrange 
these  nonzero  roots  in  order  of  decreasing  magnitude  and  they  appear  as 


h‘1’2 . 8 


where  u.*0,  u„+...+u.  =  r.  X,>...>A  . 

0*1  SL  *  1  i, 

Since  Z  ^  >  0  we  can  write  Z  Q>0  and  rk(QE^)  = 

there  exist  orthogonal  matrices  V^,  such  that 


rk(E^)=r.  Henc 


V’Q  *V  ■ 

r1  12 


0 


Sn\. I  ...  0 
1  u2 


^nX»I 
t  u 


A  (say) 


qy'  *  V'Y  j.V 
w  1  11*  2 


QY2**Vi«2V2  qY3-ViM3V2- 


'•A>.  .v.v.v  w  v  ->  ■> 


where 


\  =  /Xh  (wla8)  “  *  Vl+1 . \  . 


6  *  r+1 . u 


Then  (5.5)  can  be  written  as 


XI  . 

|n“1A  +  n“1/2C  +  .*  1*. 


d  -B 
n 


=  0  . 


(5.7) 


Let  d^,...,d  be  the  roots  of  (5.7).  We  classify  them  into  £+1  sets,  con¬ 


taining  ,  s-r  members  respectively. 


* 

For  the  last  set  let  d  be  any  one  of  them  putting  d  =  n  t  and  substituting 


it  in  (3.19),  cancelling  the  common  factor  n~  from  the  last  p-r  rows  and  columns 


and  finally  letting  n-*-«.  It  is  equivalent  to 


(5.8) 


E1  *  '  *  E£  A£+l£+ l"T  n  B£+l£+l 


where 


£+l£+l 


C  l  wltoWitB)  a, 6  -  r+l,...,p. 


£+l£+l 


n-k 

<j1"3taw3t6) 


a. 6  -  r+1 ,. . . ,p. 


Equation  (5.8)  is  equivalent  to 


lim  A. 


£+l£+l  A,  11 
n-*»  1 


-  ~  E’E,  -  . . .-  ~  E’E  -  t  -  B 


X£  £  £  n  £+!£+!' 


-  Urn  |E-  t-  B  |  -  0 

n+® 


(5.9) 


jV 
iv  ' 


*r  -  vr w  *  *i 

-v- 


-  *  - 


;-v.v,v 


m 


where 


u 

E-1  l  -VoW  “•6-r+1 . 1 


t-r+1 


Here,  we  note  that  eq.  (5.8)  is  obtained  by  following  same  lines  as  in  Hsu  (19 4 L) 


When  the  underlying  distribution  is  multivariate  normal,  we  denote  E,  By 


and  ,  respectively.  We  have 


n  f(n)  (N)"1  d  -l 

E  B £+!,£+!  nEB£+l,£+l  * 


Since 


Um  n  Crt  '  1  <»•••> 

If*00 


we  have 


lim|EW  -  T“  B  .  [ 
1  n  £+l£+l 1 

n-*« 


|E  -til  *  0  (a.e.) 


(5.10)  has  s-r  nonzero  roots,  written  as  t  . x  .  Let 

n-i  s 


x .  =  nd.  i  •  r+1, . . . ,s. 
x  i 


Then,  by  Lemma  2.3,  for  i=r+l,...,s  we  have 


-*■  x^  as  n  -*■ 00 . 


So  when  n-»-«° 


-2  In  T  -  -2  l  lnCl+d  ) 
i-rfl 


l  InCl  +  ^)n 

i-rfl  n 


*  l  t  ■  trE^-x!  w  \ 
i.r+l  1  A(p-r) (u-r) 


(5.10) 


(5.11) 


MW  tT%  a1 
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Similarly,  we  can  show  that  the  asymptotic  null  distribution  of  -2  In  is 
2 

x(p-r)  (u-r-1)  *  _n 

O.  _ 

We  know  T,  *  t(l+*  (.1+0 )  »  s  “  min(u-l.p)  and  l  >...>£->  6  are  the 

4  r+x  s  i  s 

positive  roots  of 


0  -  | Sf CHqM)  -  £S| 

-  |Y*  tly-CT"1!)  ((T'1!)  '  (t"1!)  r1®’1)  '  ]Y  -  IYJYj | 


where 


the  same  density  as  in  the  proof  of  (5.11). 


Since 


rktl^d”1!)  ((t"1!)  ’  (l”1!)  )~1(T~11)  '  ]  =  u-1 


and  it  is  an  idempotent  matrix,  there  exists  an  orthogonal  matrix  I*  u  such  that 


r [i  -(t"1i)((t‘1i),(t“1))"1(.t“1i),1  *  tV1 

u  C  C  C  C  0  0 


X  C 

Taking  orthogonal  transformation 


f 

fx.l 

r  o 

r  Y_] 

,  1 

1 

h* 

ro 

Y  * 

Y2 

Y, 

0  I 

Y0 

t  3j 

n-uj 

.  3; 

i.e. 


y1  =  r'x  y2  -  y2  y3 


=  Y. 


Partition  X  as 


X 


^X^|  u-1 


1  * 


Then  l>. .  .>£_  >  0  are  positive  roots  of 
1—  —  s 


XiXl-£Y3Y3l  *  ° 


and  the  density  of  yl  is  the  marginal  density  of 

3 


2  h[tr 


Y’Y  J 
3  3 


where 


Y2*“  W 


Since  rkCr’E^)  =  r>0.  Similar  to  (5.11),  there  exist  orthogonal  matrices  V^,V 


such  that 


v^(Qr,=1)v2 


/nX7  I  *  •  *  0 

1  u, 

*  1 

i  « 

0  »  •  •  /nA .  I 
*  u* 


►i  (*<»’ 


(say) 


where  X  Uj+...+u^»r  and  l  ^  «  Q'Q  Q  >  0.  Taking  transformations 


QX*  -  V’X*V2  QY'*  -  ViW2V2  QY3  *  V[W3V2 


and  denoting 


X*  -  wx  + 


-«>J  lo 


we  observe  that  Jl,>. . .>i,  > 0  are  the  roots  of 
1 —  —  s 


I  (w|X) +ACL))'(w|1)  +  A(1))  - 


iw'w3| 


and  the  density  of 


is  the  marginal  density  of 


6.  ASYMPTOTIC  DISTRIBUTIONS  OF  TEST 
STATISTICS  WHEN  THE  OBSERVATIONS 
ARE  INDEPENDENT 


Consider  the  model 


X  -  A6  +  E 


where  A  and  0  are  as  defined  in  C2.1),  and 

( E  * 


E 

n*p 


(1) 


E  ' 

!(n) 


But  we  assume  that  E^,..,,E^  are  distributed  independently  as 


-1/2  i 

|l|  h(trE  E(i)Ej±)) 


(6.1) 


with  characteristic  function  <t>(tr£TT').  We  discuss  the  asymptotic  distribution 

•* 

,-l 


of  the  last  s-r  nonzero  eigenvalues  of  S^(0 ,M) (Se (I,S) }  and  S^(0,M){Se(I,S)} 
Following  the  same  lines  as  in  Section  5,  we  have 


liffilE  -  T  t  bm im'  -  0 

n-«o 


(6.2) 


where 


u 


E  -  (  l  wltawlttj>  a>6  *  rfl,...,p. 


t-i+1 


n-k 


Bl+U+1  ■  CJlu3tcw3t6)  °'8  =  r+1 . . 


(6.3) 

(6.4) 


W1  “  (wla8)  "  Yl*  "  A 


U*p 


V2[Cln  °)rCClk  °^rA(X“A0)Q,]Vi 


*L‘ 


(n-k)xp - 


v2[(°  In_k)r^x-A0)Q,]V[ 


(6.5) 


and  Vi>V2  are  two  orthogonal  matrices,  Q'Q  =  > 0. 


z  =  :  =  (X-Ae)Q' 

nXP  7» 

~  (n) 


(6.6) 


From  (6.1)  we  observe  that  Z^,...,Z^  are  distributed  independently  as  h(trZZ’) 


Now,  let 


rc  -  <w-  *  <W- 

k»k 


By  simple  computations  we  have 


Jl(  JirG.lirMi)Z£;L  *  *  *  *  Jl^1Ci^irCiirAili)Z£n) 


Wl"  V2 


^l(JirCuirMi)ZJa,*'jl^  1^ircuirAili)Zin 


fz  .  .  . 


:n  ;  in 

-  V2 . v*  (say) 

•  •  •  •  • 

*  * 

Z  ....  Z 
L  ui  im 

and  for  any  a  (l,...,u),  b  (l,...,n)  we  have 


(6.7) 


Var  Z 


T“Ij'Arc.iiwv 


(6.9) 


1-li-l 


'  <w  V 

“  [-2*'(0)1  +  j.rc«irAtircbir*«) 


■  [-2*'(0)](jir^(jr^1)  +Jjrwra>J<jrAUr4tj)) 

=  -2*'(0) . 


(6.10) 


In  the  following  we  prove  that  the  different  Z  ,  '  s  are  uncorrelated.  For 

ab 

any  c(^a)  e  (1,...  ,u)  ,  d(^)  e  (1,...  ,n) 


E  ZabZcd  -  E(J<jirCairAU>Z«b)(J‘1SirCCirAU>Ztd) 

't,LlI<t<rCair“l)  (1^ircCirAmi)E(Z«bZmd)  1 
When  l  4  m:E^Z{.hZm^=  0  and  so  E^ZabZcd^  =  °*  When  1  ~  m»  b  4  d  and  E(Z„uZmJ)  =0 


£b  md 


So  E(zatjZc(p  *  0«  When  l  =  m,  b  =  d,  we  have 

E<bZ‘d>  -  <-2*,<0))XircAircc3(jEMir“J>- 


Also, 


JjA!lirAtj  "  0  1  *i 


jjL  -  fc«  jircairccl  -  0  ‘’J- 


So,  E(Z  .  Z  ,)  -  0.  Now,  let  Z  . 

ab  cd  ab 


(0)  U^.  By  the  central  limit  theorem 


we  get  that  U  ^-’s  are  mutually  asymptotically  independent  normal  variates  with 
zero  mean  and  unit  standard  deviation. 
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Similarly, 


W  =  V 
3  2 


n  n 

£  P  Z<  1  *  •  *  1  W-1Z-i, 


^  Aik+1  il  *  *  *  i^1‘Ai,k+l  ip 


n 


l  l  Wil  *  *  ‘  l  r 


i-1 


Vi 


Ain  ip  1 


Tc+1,1 


.  Z 


k+l,p 


•  • 


"nl 


.  z 


vi 


(say) 


np  ' 


and  for  any  i,k  e  (k+1, . . . ,n)  j ,1  e  (1, .. . ,p) 


E  ZtJ  -  ° 


E  zuzk* 


-2*'(0)  1-k,  j-t 

0  otherwise. 


Let 


Z±.  =  /-2*'C0)  V  . 


The  vij's  are  mutually  asymptotically  independent  normal  variates  with  zero  mean 


and  unit  standard  deviations. 
Let 


u 


E  “  (I  UtaUtg)  a, 6  -  r+l,...,p. 


t-r+1 


n-k 


B£+U+l  *  VtaVt8)  a,S  “  r+1 . . 


Then 


liffl  lE-nWll  0 

n-**> 


is  equivalent  to 


n-x» 


lim  n  BJH-U+1=  I* 
n-x=° 


>,  as  n -*■*>,  Tr+^>***»TS  are  the  eigenvalues  of  the  central  Wishrat  matrix  with 


(u-r)  degree  of  freedom. 


Then  following  the  same  lines  we  can  discuss  the  asymptotic  distribution  of 


the  last  s-r  nonzero  eigenvalues  of  S-C0 ,M){S  (I,S) }_1. 

r  e 


>v<  ' .  t- 
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